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The Fourier transform of the Hadamard transform: 
Multifractals, Sequences and Quantum Chaos 

N. Meenakshisundaram and Arul Lakshminarayan 



Abstract — We introduce a class of functions that limit to multifrac- 
tal measures and wliicli arise wlien one takes the Fourier transform of 
the Hadamard transform. This introduces generalizations of the Fourier 
transform of the well-studied and ubiquitous Thue-Morse sequence, and 
introduces also generalizations to other intriguing sequences. We show 
their relevance to quantum chaos, by displaying quantum eigenfunctions 
of the quantum bakers map that are approximated well by such measures, 
thereby extending our recent work where we pointed to the existence of 
"Thue-Morse" states. 

Keywords — Fourier transform. Hadamard transform. Multifractals. Se- 
quences. Quantum chaos. Eigenfunctions. 



I. Introduction 

THE Fourier and Hadamard transforms are standard tools, 
widely used in science and signal processing 1 1 1. The rel- 
ative importance of the two transforms may be judged to be a 
factor of thirty in favour of the Fourier transform if one were 
to go by a "google" search which returned over five million 
webpages for this transform. Both these transforms can be im- 
plemented with fast algorithms that reduce their implementa- 
tion on N data points from N'^ to N log(A^) operations. The 
fast Fourier transform (FFT) and the fast Hadamard transform 
essentially rely on the factoring of the transform into opera- 
tors acting on product vector spaces. The Hadamard transform 
though is a real transform which only adds or subtracts the data 
and is therefore widely used in digital signal processing. The 
Fourier transform conjugate spaces are familiar ones ("time- 
frequency", "position-momentum") etc., while the correspond- 
ing Hadamard transforms are not so well understood. Neverthe- 
less the Hadamard transform has also received great attention in 
the recent past due to its uses in quantum computing, with the 
Hadamard gate being a central construct |2 1. 

We define the Fourier transform (FT) on N sites as 



^=exp(-2TTi(m + a)(n + a)/N) (1) 



(Gn) 



with < a < 1 /2 being a phase parameterizing the transform, 
and < m,n < N — 1. Gn as a matrix is an unitary one. The 
Hadamard transform that we use maybe written in several ways, 
firstly as an tensor or fCronecker product form, secondly via a 
recursion and finally via their matrix elements. In all of this and 
what follows in this paper we assume N to be a power of 2, i.e., 
N = 2^, for some integer K. If 
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Equivalently 
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and Hoa = 1 . Also in terms of matrix elements 
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where a ■ b 



bi and m 

X^ili ^i2*~^, that is a and b are vectors whose entries are the 
binary expansions of the matrix positions (m, n). Note that Hn 
is such that Hfj — I, while G% — I, where / is the identity, 
therefore the spectrum of both these transforms are highly de- 
generate (±1 for Hn, ±1, ±J for Gn)- 

Also notice that we can enumerate the columns of Hn (or 
rows, as it is a symmetric matrix) as outer products of 



Wo 



V2 \ 1 
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If n — aK-iaK-2 • • • oq is its binary representation (0 < n < 
2^ — 1), the rfi^ column Vn is the outer product 



(7) 



While the first column Vb is simply an uniform string of Is, 
the last one Vn-i is the A'-th generation of the Thue-Morse 
sequence. 

The Thue-Morse sequence is an example of an "automatic 
sequence" ('31 of two alphabets say A and B. Given say A, 
the rule is to replace it by AB and given B replace it by BA. 
Thus starting with the seed A we get {A AB —> ABBA 
ABBABAAB ■ ■ ■}. The K-th generation consists of a 
string or word of length 2^^ which is cube free, that is no block 
(any finite string consisting of the two alphabets) repeats thrice 
consecutively. This sequence occurs in numerous contexts 0, 
combinatorics on words, number theory, group theory, dynam- 
ical systems, to name a few, and is considered to be marginal 
between a quasiperiodic sequence and a chaotic one. The de- 
terministic disorder of this sequence is relevant to models of 
quasicrystals |4|, mesoscopic disordered system {5j, and as we 
established recently to quantum chaos Q. The column Vn-i is 
got from the K-th generation of the Thue-Morse sequence via 
the identification A = 1 and B = —1, apart from the factor 

1/Vn- 

Our recent work indicated that the Fourier transform of the 
Thue-Morse sequence, along with the sequence itself was an 
excellent ansatz for a class of eigenstates in the quantum baker 
map, which we called the "Thue-Morse states". The classi- 
cal baker mao is a oaradismatic and simole model of comolete 
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Hamiltonian chaos. Its quantization is then of considerable in- 
terest in the study of quantum chaos, therefore the emergence of 
this ansatz provides an interesting way to think of the determin- 
istic structural disorder of quantum chaotic eigenfunctions, at 
least in this model system. We also found that the Fourier trans- 
form of some other columns of the Hadamard transform played 
a crucial role in describing other states. The Fourier transform of 
the Thue-Morse sequence 1 7 1, 1 9 1, or some of the other columns 
of i/jv not simple functions though, they could be multifrac- 
tals 1 8 1 . Of course the Fourier transform of the first column Vq is 
just a localized delta peak (which maybe broadened for nonzero 
a); thus we expect that the Fourier transform of the Hadamard 
matrix will result in a mixture of functions or measures with a 
range of complexity. This is our primary motivation for study- 
ing the product GnHj^, the Fourier transform of the Hadamard 
transform. 

II. The matrix elements of GnHn 

The matrix elements of GnH^ are evaluated economically 
as a product of K trigonometric terms. Using the matrix repre- 
sentations of Gat and we get 

1 



1=0 



(8) 

where Z = ^^^g 2^ bj, and n — J^j^o 2% j . Thus performing 
the independent sums over the bj, and after some simphfica- 
tions, we get 
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Thus we are led to the study of the following class of func- 
tions which are power spectra's of the columns of the Hadamard 



matrix: |(GArK)fcP = {{GnH, 



Njkn 



fn{k) 



K-1 



J=0 



(10) 



We view these as a function of k for a fixed n — 
aK-iaK-2 ■ ■ ■ do- They satisfy the normalization, that follows 
from the unitarity of GnHn, 



N-l 



(11) 



fe=0 



and we will in fact treat /„ (fc) as a probability measure. We are 
interested in the limit N ^ oo {or: K ^ oo). If there exist 
sequences of n that lead to limiting distributions we are espe- 
cially interested in these. In the following we use the notation 
that {s)m is an m-fold repetition of the binary string s. If n is 
of this form we also denote fn{k) as /(^^(fc). For instance the 
case when n = N — 1 = (I)a' leads to the power spectrum 
of the Thue-Morse seauence that is well-known to limit to an 



multifractal measure fT^. In this case 



K-1 



/(i)(fc)= n™' 



(12) 



This has been particularly studied when a = and found to 
limit to a multifractal with a correlation dimension D2 = 0.64 

III. The participation ratio or the correlation 

DIMENSION 

To probe the limiting functions, if they exist, for multifractal- 
ity we test for the scaling relation 



JV-l 



D2 



(13) 



fc=0 



The left hand side of this is also interpreted as the inverse par- 
ticipation ratio, its inverse, Pn, being the effective spread of the 
power spectrum, an estimate of the number of "frequencies" 
(fc) that participate in it. If £'2=0, the frequencies are local- 
ized (Bragg peaks of crystallography), if D2=l, it is a situation 
expected of power spectra's of random sequences. In the inter- 
mediate range are multifractals, with structures at many scales. 

Consider the class of functions that result as K tends to 00 
along even numbers, and n — {01)k/2- Equivalently n — 
{N — l)/3, and the functions are /(oi)(fc)- We also simulta- 
neously consider the closely related functions /(io)(fc)- Both 
these tend to multifractal measures as iiT cx3 with D2 ~ 0.57. 
The principal peaks of /(io)(fc) are at 1/5,4/5, while those of 
/(oi)(fc) are at 2/5, 3/5. Taken together these peaks constitute 
a period-4 orbit of the doubling map x ^ 2x (modi). The 
peaks of the Fourier transform of the Thue-Morse sequence, or 
of f(i){k) are at the period-2 orbits of the doubling map, i.e. 
at 1/3,2/3. We see these and a few other such functions in 
Fig. (1). The scaling of the participation ratio of these mea- 
sures and the corresponding correlation dimension are shown in 
Fig. (2), which shows that indeed these measures are multifrac- 
tals. In the case of strings of the form (001) for instance, K is 
taken to be multiples of 3 and so on. One interesting observa- 
tion from this figure is that it appears that the more 1 there are 
in the string s, the more is the dimension D2, so that the power- 
spectrum of the Thue-Morse sequence may have the maximum 
possible D2 value in this class of multifractals. Of course the 
string {i))K is not a multifractal at all, and D2 — Q va this case. 

For a given N, the participation ratio P„ of the Fourier trans- 
form for the various columns n, has a range of values that indi- 
cates the localization in the conjugate basis. We show in Fig. (3) 
the participation ratios for the case N = 1024 and a — 1/2. We 
see here the intricate way in which the columns of the Hadamard 
matrix are arranged. The largest participation ratio occurs for 
the last column of the Hadamard matrix which corresponds to 
the Thue-Morse sequence. We find a similar behaviour for other 
values of a, as well as other measures of localization such as 
the entropy. We include this parameter as for the application 
we have in mind, namely the quantum bakers map, a = 1/2 is 
oertinent. 
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Fig. 1 

THE FRACTAL MEASURES /{I), (/{oi) , /(lO) ). /(lllO). (/(OOl). 
/(OlO): /(lOO)) ARE PLOTTED CLOCKWISE. A" = 10 IN THE FIRST CASE 
AND = 12 IN THE OTHERS, WHILE a = UNIFORMLY. WE HAVE 
GROUPED SOME FUNCTIONS AND SHOWN THEM IN DIFFERENT STYLES IN 
THE FIGURE. 



Fig. 2 

The SCALING OF THE PARTICIPATION RATIO P„ WITH TV. logj (Pn) IS 
PLOTTED AGAINST logj (N) = K, AND n = (s)m , WHERE THE STRING S 
IS GIVEN IN THE KEY TO THE FIGURE, AND m IS A WHOLE NUMBER SUCH 
THAT K/m IS THE LENGTH OF THE STRING S. THE SLOPES OF THE LINES 
ARE ALSO GIVEN IN THE KEY. 



We now briefly discuss the sequence whose power spectrum 
is /(oi)(fc)- Recall that the corresponding sequence for /(i)(fc) 
was the Thue-Morse sequence {1, —1, —1, 1, —1, 1, 1, —1, • • •}. 

EK-l 
^ ,„ „. „„„ ^ i=o where 

again the a; are bits of the binary expansion of n. Similarly the 

n-th term of the sequence whose power spectrum is /(oi) (fc) is 

i„ = (-1)^-0,2.... (14) 

The first few terms of this sequence are {1,— 1,1,— 1,— 1,1, 
— 1, 1, 1, —1, 1, —1, —1, 1, —1, 1, —1, 1, —1, 1, . . .}. To write an 
concatenation rule we use the fact that this is formed by repeated 
outer product of (1, —1, 1, —1)^ and get 



S{k + 1) ^ S{k)S{k)S{k)S{k) (15) 

Where S{k) is the fc-th generation of the sequence, with S{0) — 
1, and S{k) is the complementary set where 1 and —1 are inter- 
changed. While we found similar rules and sequences elsewhere 
jlOl, we did not find this exact one. It also appears that the infla- 
tion rules A — > ABAB, B BAB A produces this sequence. 

IV. Connections to quantum chaos 

So far we have introduced the measure and discussed some 
of their mathematical properties. Here we make explicit their 
relevance to quantum chaos. The classical baker's map |11|, 
T, is the area preserving transformation of the unit square 
[0, 1) X [0, 1) onto itself, which takes a phase space point {q,p) 
to iq\p') where {q' = 2q, p' = p/2) if < g < 1/2 and 
{q' = 2q - 1, p' = {p + l)/2) if 1/2 < q < 1. The stretching 
along the horizontal q direction by a factor of two is compen- 
sated exactly by a compression in the vertical p direction. The 
reoeated action of T on the sauare leaves the ohase soace mixed. 




Fig. 3 

The participation ratio P„ vs n for N = 1024 and o = 1/2. 



this is well known to be a fully chaotic system that in a mathe- 
matically precise sense is as random as a coin toss. The area- 
preserving property makes this map a model of chaotic two- 
degree of freedom Hamiltonian systems, and the Lyapunov ex- 
ponent is log(2) per iteration. 

The baker's map was quantized by Balazs and Voros 1 12], 
while Saraceno 1 13 1 imposed anti -periodic boundary conditions, 
and this leads to the quantum baker's map, in the position repre- 
sentation, that we use in this Letter: 

C'r <><» 

where (C?Ar)„i„ = {PmWn) = exp[-27ri(m + l/2)(n + 
1/2) /N]/\/N. TheHilbert soace is finite dimensional, the di- 
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mensionality N being the scaled inverse Planck constant {N = 
1/h), where we have used that the phase-space area is unity. The 
position and momentum states are denoted as |g„) and \pm), 
where m,n = 0, ■ ■ ■ , N — 1 and the transformation function 
between these bases is the finite Fourier transform Gjv given 
above. 

The choice of anti-periodic boundary conditions fully pre- 
serves parity symmetry, here called R, which is such that 
R\qn) = |(7Ar_„_i). Time-reversal symmetry is also present and 
implies in the context of the quantum baker's map that an over- 
all phase can be chosen such that the momentum and position 
representations are complex conjugates: GN<f> = 4>* , if </> is an 
eigenstate in the position basis. B is an unitary matrix, whose 
repeated application is the quantum version of the full left-shift 
of classical chaos. There is a semiclassical trace formula, which, 
based on the unstable periodic orbits, approximates eigenvalues 

ini. 

The nature of quantum chaotic eigenfunctions is intriguing as 
they display a bewildering variety of patterns, that can some- 
times be partially attributed to classical periodic orbits. This 
phenomenon called scarring 1 15| is apparently in conflict with 
another observed behaviour namely their similarity to random 
matrix eigenfunctions 1161 . The quantum bakers map affords 
us an opportunity to study such states in a simple setting, how- 
ever no known analytical formulae exist. We had proposed a 
few ansatz for a variety of states based upon the Thue-Morse 
sequence and its Fourier transform f6l. These could sometimes 
reproduce states to more than 99%. In this work we point to 
other measures such as /(oi) that also play a role in the spec- 
trum of the quantum bakers map. 

In Fig. (4) we see two examples of states of the quantum bak- 
ers map in the position basis for N = 1024, along with their 
Hadamard transforms. The Hadamard transforms are particu- 
larly simple, the Thue-Morse state, iptm, being the first |6|. We 
show it here for comparison with another state whose Hadamard 
transform has distinct peaks at around iV/3 and 2iV/3 implying 
that it is likely that the measures /(oi) /(lo) discussed above 
are relevant for these states. There are usually more than one 
such state, we show here a particularly "clean" state, in terms of 
its Hadamard transform, and for the purpose of this paper call it 

'0(01)- 

Indeed we can construct the ansatz 

0= (7l+7rG'-')^(iV-l)/3 + (72 + 72G~')l^2(Ar-l)/3 (17) 

where N are powers of 4. This includes both the columns at 
(A'^ — l)/3 and 2{N — l)/3, along with their Fourier trans- 
forms. The latter are included due to the presence of time- 
reversal symmetry. In fact they dominate the state and is the 
motivation for our study. This ansatz has two complex constants 
7i and 72 which we determine numerically so that its overlap 
with the actual state ^/'(oi) is the maximum. For instance in 
the case N — 64, we were able to find 7 such that the over- 
lap |(</'|'0(oi))P ~ 0.75. A comparison of the spectral mea- 
sures for the Thue-Morse sequence, as well as /(qi) and /(ip) 
in Fig. (1) with the actual eigenfunctions in Fig. (4) show the 
similarity between them. Notice that the actual wavefunction is 
more closely related to the Fourier transforms of a linear com- 
bination of the columns of the Hadamard matrix. The oeaks of 
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Fig. 4 

The Thue-Morse state V'tm (top left) and a state called in the 

text l/'(01) (bottom left), along with their corresponding 
Hadamard transforms (right). In this case N = 1024 and the 
peaks of the hadamard transform of l/'{01) occur at (tv — l)/3 

AND 2{N ~ l)/3. 



this wavefunction can be identified with the period-4 orbit start- 
ing from [q = 1/5, p — 4/5), and has been noted earlier the q 
part of this corresponds to the peaks in the measures /(oi) and 

/(lO)- 

We have concentrated on convergent measures as n — (s)™ 
tends to infinity and s is a bit string repeated m times. We can 
have other sequences of n that lead to convergent measures, in 
particular those whose binary expansions ends in a string of 
1 or 0. In the former case we have already used such mea- 
sures to describe a family of states of the bakers map that are 
scarred by period-2, period-4 and associated homoclinic orbits 
||6l. As spectral measures they are different depending on the 
finite string that precedes the infinite string of 1 ; however their 
correlation dimension D2 seems to be the same as that of the 
spectrum of the Thue-Morse sequence, namely 0.64. As we 
noted earlier this seems to be consistent with our view that more 
the 1 , more the dimension D2 ■ If the string ai is itself not peri- 
odic, say when n is tending to infinity such that tends to 
an irrational number, there does not seem to be any convergent 
measures. 

In summary, we have introduced a class of simple functions 
that limit to multifractals, have interesting connections to se- 
quences and to a simple model of quantum chaos. We have done 
this by combining two well-known, well-studied and standard 
transforms, namely the Fourier and the Hadamard. 
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